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Abstract 

We make a numerical study of the classical solutions of the combined sys- 
tem consisting of the Georgi-Glashow model and the SO(3) gauged Skyrme 
model. Both monopole-Skyrmion and dyon-Skyrmion solutions are found. A 
new bifurcation is shown to occur in the gauged Skyrmion solution sector. 
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1 Introduction 



There are two 3+1 dimensional SO (3) gauge field models which support static soliton 
solutions. One is the Georgi-Glashow (GG) model which supports the well known 
monopole and the other is the SU(2) L+R , or vector, gauged Skyrme model [Q, [| 
which also supports £0(3) gauged Skyrmions. In addition to the monopole, the GG 
model supports also dyon solutions 0] which in addition to the magnetic charge carry 
an electric charge as well. The topological stability of the monopole comes from the 
magnetic charge, which is descended from the second Chern-Pontryagin charge, while 
the topological charge of the gauged Skyrmion is the degree of the map. 

Combining these two models, we have a new system whose topological charge is 
a sum of the respective charges, and it can reasonably be expected that this system 
also supports static finite energy solitons. Note that in this case the local SO (3) 
symmetry is broken down to U(l) via the Higgs mechanism, in contrast to the 5*0(3) 
gauged Skyrme model on its own, in which case the local £0(3) symmetry is not 
broken at all and 3 massless gauge bosons survive. In this preliminary investigation, 
this is precisely what we have done. Using numerical methods, we verify that such 
solutions exist. Moreover, we have sought and found both monopole-Skyrmion and 
dyon-Skyrmion solutions, and studied some of their properties. The combined system 
supports solutions also with zero monopole charge, unit Baryon charge, as well as with 
unit monopole charge, zero Baryon charge. 

Even though this is a self contained numerical study of the classical solutions 
alluded to above, it is in order to put it into context both in the background of 
previous work involving the gauging of the Skyrme model || , and, from the viewpoint 
of its potential physical relevance. 

The Skyrme model was gauged by Witten in Ref. ||, and others e. g. in Refs. 
0. These works were carried out in the context of current algebra results, and were 
not concerned with the solitonic aspects of the gauged Skyrmion. That was done 
subsequently by many authors |TJj where gauged Skyrme solitons were studied with 
the aim of explaining the low energy properties of Hadrons. Also in the context of 
electro weak theory, which can be regarded as a gauged Skyrme model in the limit of 
very high Higgs mass, Rubakov || and Eilam et al. [|10| considered the static classical 
solutions of the SUl gauged Skyrme model. In all these cases, there is no topological 
lower bound and the classical solutions are metastable, but for certain values of the 
parameters in one of these models || [K| a stable branch of solitons appears as a 
result of catastrophic behaviour. The advantage of the gauging used in Refs. || |TU 



is that the 4-divergence of the topological current does not vanish but equals the local 
chiral anomaly [|ll| |l2j, which can present itself as a mechanism for Baryon number 
violation as explained in Ref. 0. 

In the context of Baryon number violation, there is an older mechanism suggested 
by Rubakov [13j and by Callan |14| where monopole- (left-handed massless) Fermion 
interactions lead to Fermion number nonconservation. The mechanism involves the 
fluctuations of the electric field, in the presence of the magnetic field of the monopole, 
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giving rise to nonzero chiral anomaly and hence Fermion number violation. This 
was shown for the case of massless (left-handed) Fermions, by scattering with the 
monopole, which describes a high energy process. The approximation techniques 



employed |14], [13[] are neither perturbation theoretic nor semiclassical. To describe 



a low energy process such as a decay, it would be more appropriate to deal with a 
process that is susceptible to semiclassical analysis. To this end, Callan and Witten 
fl5f replaced the massless Fermions by the Skyrme soliton ||, interacting with the 
(Abelian) magnetic field of the monopole. While they |H| did not seek to demostrate 
the existence of a U(l) gauged Skyrmion, this is implicit in their work and has 
recently been verified numerically |T7|]. In the background of this it is hoped that 
the present work, which sets out to find the monopole-Skyrmion and dyon-Skyrmion 
solutions, would be of concrete usefulness to a semiclassical method of describing 
Baryon number decay. In particular the dyon-Skyrmion excites a nonzero classical 
quantity for the chiral anomaly, which can lead to chirality breaking as pointed out 
long ago by Marciano and Pagels [M . 

In Section 2 we present the model and give the topological lower bounds on the 
static energy. In Section 3 we give the static spherically symmetric fields and the 
field equations in the static limit. Sections 4 and 5 deal, respectively, with the results 
of the numerical analysis of the Ao = and Ao ^ cases. Section 5 in particular, 
includes an in depth analysis of the Julia-Zee dyon J|. We summarise and discuss 
our results in Section 6. 



2 The Model 

The model under consideration is the combination of the Georgi-Glashow (GG) model 
and of the 5*0(3) gauged 0(4) (Skyrme) model studied previously in Refs. [§, ||. We 
state the Lagrangian of each of these models separately, defined in four dimensional 
Minkowski space, each being normalised properly so that the value of the energy of 
the static soliton in each case lies above its own topological lower bound. The static 
solutions in question satisfy the Euler- Lagrange equations of the static energy density 
functional, which is the static hamiltonian in the temporal gauge. In the GG case, 
this is the 'tHooft- Polyakov |l| monopole, while in the latter case it is the soliton 
studied in Refs. |2|, |3|. 

The GG model is described by 

Ccc = -\^\ 2 + \xi\D^ - \\i(v 2 ~ l$1 2 ) 2 , (1) 

F« v = d II A:-d u A a ^+e a ^A^Al D^ a = d^ a + e a ^A^\ (2) 

The late Greek indices /z, z/, .. label the Minkowski space vectors, while the early Greek 
indices .. = 1,2,3 label the elements of the algebra of the gauge group S'0(3). 
The Latin letters a, b, .. = 1, 2, 3, 4, so that a = (a, 4) are reserved for the 0(4) Skyrme 
model. In eq. ([I]) the constant rj is the vacuum expectation value (VEV) of the Higgs 
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field and like the latter has the inverse dimension of a length. The constants Ao, Ai 
and A 2 are all dimensionless. 

The SO (3) gauged Skyrme model is described by 

£o ( 4) = -\4\FP 2 + ^?IVI 2 - l4\D[^ a D^ b \ 2 (3) 



D^ a = d^ a + e a ^A^\ D^(f) = (4) 

The constants k,q and K2 are dimensionless and the constant K\ has dimension of a 
inverse length. The reason we keep all the coupling constants arbitrary in eqs. (|l|) 
and (|3|) will appear soon. 

When we consider the static Hamiltonians corresponding to the Lagrangians above 
in the temporal gauge, i. e. A = 0, we can write the following topological identities: 

J dr H GG > A^XlXlM , (5) 

n GG = \\ 4 \F^ + l -\\\D^\ 2 + - |$"| 2 ) 2 , (6) 

where the integer M, representing the index of the mapping is the monopole 

topological charge. Similarly || 

/ dr H u) > 127r 2 Ki«l— p 1 T ( 7 ) 

WoW = \4\F« k ? + + l -4\D [3 ct> a D k] 4> b \ 2 (8) 

where the integer T, representing the index of the mapping (f) a (x), is the Skyrmion 
topological charge. In the Skyrme description of hadrons, T is identified with the 
baryon number. 

In the following will consider also the equations resulting from the superposition 
of the two lagrangians eqs. ([l]) and 

£m = I~-gg + £0(4) • (9) 

The finite energy configurations of this mixed lagrangian are characterized by the 
couple of integers M, T. Classical solutions corresponding to the two different topo- 
logical excitations can then be constructed, they correspond to the configuration with 
minimal energy in a class M, T. 

In order to normalize the fields conventionally, we have to choose 

\ 2 = -cos(9) , K 2 = -sm(9) , \* = 1 (10) 

e e 

where e denotes the gauge coupling constant. With the choice 9 = tt/A, the topological 
inequality relating TC m to the class of solutions of indexes M, T reads 



where 

a = 4 , t = -A , * = e 2 4- (12) 

In eq. fll2|) , A is defined for later convenience. Note that the topological lower bound 
eq. ( |iT|) can be refined by an optimal value of the mixing angle 9, depending on the 
parameters Ai,A2, Ki,K2- To acheive this it is necessary to solve a complicated non 
linear equation, which we shall not pursue here. 

3 Static spherically symmetric equations 

The classical equations corresponding to eqs. ([!]), ([|) and (JS)) are in general intractable. 
We will restrict our search of solutions to the static and spherically symmetric case. 
If we choose to employ the temporal gauge in the static limit, the Euler- Lagrange 
equations will reduce to the variational equations arising from the static Hamiltonians 
pertaining to the Lagrangians eqs. ([I]) and (|3|). The latter would be bounded from 
below by the monopole charge and the Baryon number densities, respectively. Hence 
the solutions to the classical equations of each of these static Hamiltionians, sepa- 
rately, can describe the 'tHooft-Polyakov monopole [jXJ] and the soliton of the 5*0(3) 
gauged Skyrme model 0, [!|. The Euler-Lagrange equations of the Hamiltonian of the 
combined static system, i. e. GG-Skyrme, in the temporal gauge also supports soliton 
solutions since the Hamiltonian is again bounded from below by the two topological 
charges eq. ( [LT]) • This is one of the problems studied in the present work yielding the 
monopole-Skyrmion solitons. 

If instead of employing the temporal gauge we proceed like Julia and Zee and 
solve the Euler-Lagrange equations pertaining to the Lagrangian eq. @ defined on 
Minkowski space in the static limit, the resulting solutions of the GG-Skyrme system 
describe the dyon-Skyrmion. This is the other problem studied in this work. As in 
the case of the dyon || on its own, we shall restrict to the spherically symmetric 
solutions only []. 

The spherically symmetric Ansatz employed is 

A" = SiapX^ , Aq = — — x a (13) 

$ Q = rj h(r) x a (14) 

<j) a = sin f(r)x a , 4 = cos/(r). (15) 

Notice that the functions a(r), h(r), g(r) and /(r) are dimensionless. We find it 
useful to introduce a dimensionless radial variable 

x = M w r , M w = erj (16) 



1 In this case the classical equations simplify sufficiently to become tractable. To our knowledge 
the only dyon solutions known are the spherically symmetric Julia-Zee [0 dyons. 
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Substituting the Ansatze eqs. flT3|)-(|T5D into the static limit of the Lagrangian 
eq. @, leads to the following one dimensional (radial) Lagrangian density L m , defined 
by 

j L m dx = Jc m dr = E 1 ~ E 2 (17) 



with 



An ~ An f 

E p = — r/Ep = — V / dx8 p > P = !> 2 
e e J 



£ 2 = (aT4 (a2 - 1)2 



2a; 2 

+ - a ; 2 (/i / ) 2 + a 2 /i 2 + ^ 2 (/ i 2 -l)' 
2 4 

+ |x 2 (/') 2 + 2a 2 sin 2 /] 



sin 2 / 

2x 2 



+ Ka 2 sin 2 /[(f) 2 + a 2 — (20) 



where we have separated the contribution E% due to the electric field and the prime 
denotes the derivative with respect to x. 

The static classical equations corresponding to the Lagrangian density C m , in the 
spherically symmetric Ansatz, turn out to be equivalent to the equations obtained 
by varying the effective one dimensional density (see eqs. (0) and (|20D) £\ — £2 
with respect to the radial functions a,g,h and /. These equations are obtained 
straightforwardy and we do not list them here. We note however, that for each 
function the corresponding variational equation can be solved trivially by setting this 
function to zero. 

It will be useful to present their asymptotic forms in the i>1 region, to facilitate 
subsequent explanations. They are, in order of the variations of a, g, h and / : 

a" =a(^-^ + h 2 -g 2 + ^sm 2 f + ...) (21) 

(x 2 gj =2ga 2 (22) 
(x^hj =h(2a 2 + \x 2 (h 2 - I)) (23) 
(x 2 fj =2a 2 sin/cos/ + o( K /0. (24) 

Following 0], we define the energy of a configuration by E = E\ + E 2 , which 
coincides with the volume integral of the static Hamiltonian obtained in the usual 
way from the gauge invariant stress tensor. The topological lower bound for E 2 follows 
immediately from eqs. (0), ([?[) and (PI). 
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4 Numerical results, case Aq = 



We first discuss the classical solutions in absence of the electric field, i. e. with 
g(x) = 0. Eq. (f22) is trivially solved and we are left with a system of three non linear 



differential equations. Only the part E 2 of the action is relevant in this case. In the 
following, we will conveniently denote the value E 2 of the solution with given M and 
T by 

E MT {\t,K). (25) 
We now describe the four cases with M < 1 and T < 1. 

4.1 CaseM = 0,T = 

This corresponds to the class of the vacuum which is not spherically symmetric. It 
has a zero energy 

£oo(A,£,/c) = 0. (26) 

4.2 Case M = 1, T = 

This case corresponds to the celebrated SU(2) magnetic monopole [|1[|. Since T = 0, 
it has f(r) = 0; as a consequence, the parameters £ and k are irrelevant for this case. 
The boundary conditions and asymptotic behaviour of the functions a, h read 

a(0) = 1 , h(0) = (27) 
a(x) ~ Ae" x , h(:r) ~ 1 - Be'^ (x -> oo) (28) 

where A, 5, F are constants. The values of the energy of the monopole solution 
were computed long ago [Tj| (our numerics fully reproduces these values); the energy 
increases monotonically with A as demonstrated in Table 1. 

In the Bogomol'nyi limit, A = 0, the energy coincides with the topological lower 
bound, i. e. (omitting the parameters £ and k) 

E 10 {X) > E 10 {0) = 1 . (29) 

The solution, the Prasad-Sommerfield monopole, is expressed in terms of elementary 
functions f[9" |. Its behaviour near the origin is given by (B^) but, for x — > oo, we have 



a(x) ~ xe x , /i(x) ~ 1 (30) 

x 



instead of (EST). 
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4.3 CaseM = 0,T=l 

In this case, there is no Higgs field and consequently no Higgs potential. Correspond- 
ingly, the topological lower bound reads 

E Q1 (\,t,K)>^^§= (31) 
2 V 1 + 9k 

independently of A. The static equations describe the gauged Skyrmion studied in 
Ref. i. 

The parameter £ can be changed by a rescaling of the radial variable x. Compar- 
ison of the energies of the gauged Skyrmion and of the monopole is demonstrated in 
Table 1 and Fig. 1. In Table 1 we list the energies for various values of A for fixed 
£ = 1 and k = 0.4, while in Fig. 1 k varies and we have fixed A = £ = 1. 

Let us now come to the detailed discussion of the solutions in the region k ~ 0.8. 
For completeness, it is useful to summarise the possible boundary conditions available 
for the gauged Skyrmion. At the origin x = the behaviour of the radial functions 
is uniquely determined by the condition of continuity of the fields at the origin: 

a(x) = 1 + A x x 2 + o(x 3 ) , f(x) = it + Fix + o(x 2 ). (32) 

In contrast, in the x ^> 1 asymptotic region, several conditions are consistent with 
the finiteness of the energy. Classical solutions of the equations have been obtained 
with the two following sets 

Type A: a « 1 - - , / « ^ ( 33 ) 

x x z 

Type B : a « 4 > / « ~ ( 34 ) 

x a x 

where a = {y/AF 2 - 3 - l)/2. 

The following results were obtained in ||. For small values of k, the solution is 
of type A, its energy increases monotonically from E = (for k = 0) and the branch 
(say branch A) stops at a critical value K — ^ A. 

~ 0.8091. For large values of k (in 
fact for k > ~ 0.69122) the solution is of type B. We call this branch B. By using 
arguments of catastrophe theory fL0 |, one can reasonably expect the occurence of a 



third branch of solutions on the interval k G [k^, Kj[], as was explained in 0. 

A third branch indeed exists. The solutions on this branch obey the condition of 
type A and therefore we refer to it as branch A. The energies of the three branches of 
solutions are depicted on Fig. 2. The branches A and A terminate at k = k c J, forming 
a cusp catastrophy. The transitions of the profile of the solutions from branch A to 
branch A is smooth. 

In contrast, when the limit k — > k c ^ is considered, the solutions of the branch A 
approach the limit of branch B in a subtle way. For instance, the value x m for which 
the function a(x) has a minimum (say a m ) tends to infinity, while a m tends to zero. 
For values of k close to Ki, the solutions of branches A and B coincide on a large 
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interval of x (typically on x G [0, 10 7 ] for k = 0.6914) and deviate from each other 
for larger values of x. In the limit k — > this interval becomes infinitely large and 
the two solutions deviate at infinity. This can clearly be seen from Fig. 3. A similar 
demonstration can be made for the function f(x), namely that near the critical point 
k c £ these functions for the two solutions on brances B and A also coincide. We do not 
display the graphs analogous to Fig. 3 in this case. The behaviour of the solutions is 
further illustrated by Figs. 4 and 5 where we plot respectively the value of F\ (defined 
in eq. ([32])) for the three branches and the value of a as a function of k. 

Fig. 2 furnishes a simple interpretation of the three solutions. To discuss it, we 
introduce the value of k where the energy of the branches A an B coincide 

(^ab ~ 0.785). On the interval k G , k c £ b ] the solution on branch A constitutes the 
absolute minimum of the energy functional E2, while the one on branch B is a local 
minimum. The solution on the branch A is a sphaleron corresponding to a saddle 
point which represents the energy barrier between the two minima. The situation 
is similar on the interval K G [Kj e ,«;J]; the absolute (resp. local) minimum energy 
configuration is then on branch B (resp. A). As k approaches the critical value k% the 
local minimum of branch B approaches the sattle point of branch A. At the critical 
k both coincide and form an inflection point. For k < k% this point is no longer an 
extremum and the solutions of branches B and A cease to exist. The global minimum 
of branch A is then the only extremum and only branch A solutions exist. The same 
scenario applies at the other critical value where the solutions of branches A and 
A stop to exist and only branch B solutions exist. 



4.4 Case M = 1, T = 1 

It is natural to call this solution the "monopole-Skyrmion". The three functions 
a, h, f are non trivial and obey the following boundary conditions at x = and as 
x — > 00, repectively, 

a(0) = 1 , h(0) = , /(0) = 7T (35) 

a(x) ~ Ae~ x , h(x) ~ 1 - B e -^ x , fix) ~ — (36) 

x 

where A, B and F are constants. In contrast to the case of the gauged Skyrmion 
solution , the finite energy condition leads to a unique asymptotic behaviour of the 
solution. 

The energy of the solution is given in Table 1 for several values of A (for £ = 1 and 
k = 0.4), indicating that the energy of the monopole-Skyrmion varies rather slightly 
with A. The corresponding lower bound inequality reads 

E 11 (X,^k) > ^ + — — r ^S= . (37) 
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4.5 General properties 

We have constructed numerically the three non trivial topological solitons above for 
numerous values of the coupling constants A, £, k and computed their energies. In 
order to give an idea of the relative magnitudes for the different classes, let us choose 
A = l,f = 1,« = 0.4, then 

E 00 = , E 10 w 1.29 , E 01 w 2.98 , E u w 3.53 . (38) 

The Bogomol'nyi limit A = is of particular interest since in that case the monopole 
saturates its topological lower bound. Choosing again £ = 1, k — 0.4, we find 

E 00 = , £ 10 = 1.0 , £ i~2.98 , E u « 3.45 (39) 

The behaviour of the solutions in the limit k — ► 0, with A, £ fixed, was carefully 
analysed. Our numerical analysis strongly supports the following formula : 

\imE m (\,£,K)=E MO (\,Z,0) , for M = 0,1 (40) 

as illustrated by Fig. 1. Indeed, in the limit k — > 0, the functions a(r),h(r) repre- 
senting the solutions of the M = T = 1 sector approach the profile of the monopole 
solution (i. e. M = 1,T = 0). At the same time, the function f(r) is more and more 
peaked at r = (in particular lim K ^ /'(0) — °°) an d tends to zero if r ^ 0. 

This result demonstrates in particular that the coupling of the Skyrmion to a 
monopole cannot stabilise the Skyrmion; the Skyrme term is necessary to guarantee 
a localized structure to the T = 1 soliton. 

The same phenomenon occurs with the branch of the gauged Skyrmion (M = 
0, T — 1) ||. The energy in this limit tends to zero, namely to the energy of the 
vacuum (M = T = 0). 

A remark should be made concerning the interpretation of the monopole-Skyrmion 
as a bound system of a monopole with magnetic charge M and a gauged Skyrmion 
with baryon number T. Consider a monopole located in a region U m centered at 
a point x m and a gauged Skyrmion located at in a region Usk centered at a point 
xsk far away from each other. Then the Skyrmion field and the corresponding gauge 
field will vanish outside the region Usk- Consequently, U m contains a pure monopole, 
consisting of a gauge field and a Higgs field. Outside U m the gauge field will vanish, 
however, the Higgs field does not vanish. Instead it will be equal to its VEV ($) va c- 
In the region £/"sk containing the Skyrmion the non- zero Higgs field is still present and 
we have to allow for interaction with the gauge field, \Dj(Q)® ac \ 2 . The Higgs vacuum 
is a constant far away from the monopole and generates masses for the gauge fields. 
Furthermore, the Higgs vacuum breaks the rotational symmetry. Consequently, the 
gauged Skyrmion solutions in the presence of a constant Higgs field will no longer pos- 
sess spherical symmetry. In addition, it might be expected that the electro-magnetic 
flux will not vanish. If we impose the condition for the Higgs vacuum that the inter- 
action with the gauge field has to vanish, Dj($)" ac = e a/37 A^($)^ ac = 0, then we will 
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find that the gauge field has to be parallel to the Higgs vacuum in isospace. This also 
breaks the spherical symmetry. 

To conclude, the interpretation of the monopole-Skyrmions as a bound state of a 
spherically symmetric monopole and a spherically symmetric gauged Skyrmion seems 
to be misleading. 



5 Numerical results, case Aq 7^ 



In order to obtain a non trivial function g(x) from eq. (|22), a non vanishing asymptotic 
value, say q, for this function has to be imposed ||. In the asymptotic region x ^> 1 
eq. ( P%D is satisfied by 

g(x)=q-^ + o(x~ 2 ) (41) 



where q, c\ are constants and q plays a major role in the construction. The equations, 
together with the finite energy condition require < q < 1, which can be seen as 
follows. In eq. (Ell) the Higgs field and the dyon field contributions, h 2 (x) — g 2 (x), 
generate asymptotically the mass term m? s = 1 — q 2 for the gauge field function 
a(x). For q > 1, m 2 a ^ becomes negative and leads to an oscillating function a(x) in 
the asymptotic region. Consequently, the term a 2 g 2 in eq. (|17]) is not integrable and 
no dyon solution exists for q > 1. 

The electric charge, as defined in M, is directly related to the constant C\ : 



Q 



4 . p oi dSi = -Q 
Attti J e 



1 

47T 



,dg_ s 
dr' 



sm 



"Cl, 



(42) 
(43) 



having used the Ansatz (13)-(|14) and eq. (f4T|) . 

Another very interesting quantity is the chiral anomaly due to the dyon-Skyrmion 
soliton whose classical solutions will be studied numerically. The anomaly equation 
for the chiral charge is 



dt 



8tt 2 



dx E,- • B,- 



f- 2 ^[g(r)(a(r)-l)}Z = fq, 

07T Z Z7T 



(44) 



having used the Ansatz ([T3|)- ([T^) and, eq. (^l|). We now discuss the solutions by 
adopting the same presentation as in the previous Section. 



5.1 Case M = 1, T = 



Here we present an in depth analysis 



The solutions are the dyons of Julia and Zee 
of this solution. The limit A = corresponds to the Prasad- Sommerfield dyon |L9 
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(PS dyon). It is worth analysing this case separately because the solution can be 
computed analytically and it provides a good check of our numerical routines. 



Case A = 



The profile of the radial functions of the PS dyon reads |y| 



a(x) = . l X (45) 
smh{cx) 



g(x) = n =^(cothcx -) (46) 

h(x) = 7= = 2 (cothca:-— ) (47) 



, = ;) (cothcx 

V 1 — q cx' 

c / , 1 
(cothco: 

\/l — cx' 

and the PS monopole is recovered for q = 0. Our parameter q is related to 7 of 
Ref. [n|, by g = tanh(7). We have chosen the arbitrary scale in the PS solution 
c = a/1 — g 2 so that the asymptotic value of the Higgs field function h(x) of the PS 
solution given above be equal to 1, since we are also studying the dyons for the A > 
case where the asymptotic value of h(x) equals 1. The charge and energy of the PS 
dyon are given by 

1 s 1 



* = 57^? ' a = (1 " 5 92) 7T^ K 1 + + o( ^ 6) (49) 

For small values of q the "magnetic" contribution to the energy, E 2 , varies slightly 
with q, accounting for the feed back of the electric charge on the classical magnetic 
energy. We would like to stress that our numerical results are in full agreement with 
these exact formulas. 

The dependence of the charge Q of the PS dyon as a function of q is represented 
in Fig. 6 (curve a). Similarly we have reported on Fig. 7 (curve a) the energy E of 
the PS dyon as a function of Q. Clearly the energy and the charge of the PS dyon 
can be arbitrarily large when q — > 1. 



Case A ^ 



For the dyon solution, the boundary conditions for the function g(x) can be read 
from eq. flU]) , and those of the functions a(x), h(x) from eqs. (|27"D and (p8|) , with the 
exception of the behaviour of the function a(x) in the x ^> 1 region, which now takes 
the form 

a(x) ~ Ae-^-^ x . (50) 

The main distinguishing feature of the A 7^ dyon vs. the PS dyon is that its electric 
charge and its classical energy are bounded for q 6 [0, 1]. 

This phenomenon appears clearly in Fig. 6 and Fig. 7 respectively, where the 
quantities Q as a function of q, and E as a function of Q, are plotted for A = 0.5. 
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More generally, it appears that the electric charge of the dyon constructed with a 
given value of the parameter q decreases when A increases. The three bullets on 
Fig. 7 represent the data given in |Q; according to our numerical results they should 
lie on line b. Our numerical results therefore slightly desagree with [Q. 

The star on line b of Fig. 7 indicates the maximal accessible charge of the dyon 
solutions for a fixed value of A. This contrasts with line a which asymptotically tends 
to infinity, in agreement with eqs. fl4gj) and fl49|) . The solutions with maximal electric 
charge and energy correspond to the case q — 1 which we discuss next. 



Case q — 1 



In the limit q = 1 in eq. ( f4T|) the equation (|21|) ceases to impose the exponential 
decay eq. (|50|) for the function a(x); we have instead 

a(x) ~ Ae~ y/ ^* for x — > oo (51) 

where ci is defined in eq. ([4lD . 

Fixing A 0, the electric charge (and similarly the classical energy) of the dyon 
cannot exceed a critical value, say Q C r(A). This quantity is plotted against A on Fig. 8 
(solid line). 



5.2 CaseM = 0,T=l 

No finite energy dyon-like solutions supporting a non- vanishing (non-Abelian) electric 
field can be found in this case. Due to the absence of the Higgs field (h = 0), eq. (f21|) 
leads to an oscillating asymptotic behaviour of a(x). The term a 2 g 2 in the energy 
eq. ( p~T|) can therefore not be integrated. 



5.3 Case M = 1, T = 1 

The boundary conditions compatible with a finite energy solution in this case are 
identical to eqs. (|35"D , ( |36D and (flTl). It is possible to construct the dyon-Skyrmion 
solutions. The dyon-Skyrmion display many features of the dyons, discussed at length 
above. These features are illustrated by Figs. 6 and 7 (dashed curves c, d and e) and 
by Fig. 8 (dashed line). In addition we illustrate the dependence of the energy on 
the Skyrme coupling constant k in Fig. 9 (solid line) for q = 0.5 and A = 0. The 
energy is an increasing function of k. In the limit of vanishing k the energy of the 
dyon-Skyrmion converges to the energy of the dyon-monopole. This can be compared 
with the behaviour of the energy of the monopole-Skyrmion, shown (for A = 1) in 
Fig. 1, where for vanishing k the energy tends to the energy of the monopole. against 
support a For our considerations leading to our conclusions in Fig. 9, we have chosen 
q = 0.5 as a typical value in the allowed range < q < 1. We expect that our results, 
summarised by the solid curve in Fig. 9, is typical for any allowed value of q, and 
also for any value of the Higgs coupling constant A, except in the important case of 
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A = and q = 1. The dyon-Skyrmion characterised by the boundary value q = 1 in 
the A = model has peculiar and interesting properties which we analyse in the next 
paragraphs. 

For A = the solutions of eqs. (^) and (|23|) are proportional to each other. 
Assuming that h — 1 at infinity, the proportionality constant is given by q, eq. (fXTl) . 
Thus for q = 1 the functions h(r) and g(r) are identical. In this special case h 2 (r) and 
g 2 {r) cancel each other in eq. fl21~|) . Consequently, eqs. ( fllD and (|4j) reduce to the field 
equation of the gauged Skyrme model (Sec. 4.3), and can be solved independently of 
h(r), g(r). 

The solutions of these equations are now given by the branch B solutions of the 
gauged Skyrme model. Once a solution for the function a(r) is found the equations for 
the functions h(r) and g(r) can be solved. Recalling that the branch B solutions exist 
for all k > K c g, we expect the existence of the dyon-Skyrmion solution for the same 
range of coupling constants k. However, not all of these solutions are finite energy 
solutions. This can be seen easily by inspecting the static Hamiltonian 8 = 81+82 
given in eqs. ([19]) and (^0|) , where the contributions from the functions h(r) and g(r) 
do not cancel. The asymptotic behaviour of these terms is dominated by a 2 (r)h 2 (r), 
a 2 (r)g 2 {r). Using the boundary conditions h(oo) = 1, g(oo) = 1 and the asymptotic 
form of the function a(r), eq. fl34"|), these terms behave like ~ A 2 (k) / x 2a ^ for large 
x, where a{n) is a function determined numerically. Thus the integration of these 
terms will give finite contributions only if a(n) > 1/2. This restricts the range of the 
coupling constant k to Ky 2 < K < 00, where Ky 2 is defined by a(/«f/ 2 ) = 1/2. For 
f = 1 we find k c { /2 = 0.7652. 

In Fig. 9 we show the dependence of the energy on the coupling constant k for 
q = 1 (dotted line) and for q = 0.5 (solid line). (For q < 1, as stated above, solutions 
exist for all values of k, the energy is a monotonic function of k, and the limit k — > 
the energy approaches the energy of the dyon solution.) For q = 1 the energy is an 
increasing function of k for large values of k only. It has a minimum at k = 1.21. As 
n approaches its critical value K,y 2 the energy becomes increasingly large and diverges 

The charge Q of the solutions is determinded by the asymptotic behavior of the 
function g(r), eq. (|43|) . Solving the equation (|^) for large x we find the following 
expressions for the charge 

lim I ci — — x _(2q_1) I for a > - , 

\ 2a - 1 J 2 

Q = \ lim (2A 2 \n(x)) for a = -, (52) 

2^00 v / 2 

lim (- J^L^-iA for a < I . 
V 2a - 1 J 2 

Thus solutions with finite charge exist only for a > 1/2, i. e. for the same range of 
the coupling constant k where finite energy solutions exist. 

In Fig. 6 we show the dependence of the charge on the parameter q for k = 0.4 
and k = 1.0. For k = 0.4 (< Ky 2 ) there is no finite charge solution for q = 1. 
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Consequently, the charge as a function of q diverges as q approaches the value 1. In 
contrast, for k — 1 (> Ky 2 ) the solution with q = 1 exists and the charge is finite for 
all values of q G [0, 1]. 

In Fig. 7 the energy as a function of the charge is shown for k = 0.4 and k = 1.0. 
For k = 0.4 (< Ky 2 ) the energy and the charge can take arbitrarily large values. In 
this case the energy is a monotonically increasing function of the charge with no end 
point. For k — 1 (> ^i/ 2 ) the energy is again a monotonically increasing function of 
the charge. However, only finite energy and charge solutions exist for this value of k. 
Thus the graph of the function E(Q) ends at the maximal value of the charge. 

6 Summary and Discussion 

We have found monopole-Skyrmion and dyon-Skyrmion solutions to an SO(3) gauged 
Higgs and 0(4) sigma (Skyrme) model, in which both scalar matter fields interact 
with the gauge field but not with each other. The Higgs field is isovector, like in the 
GG model, while the S 3 valued (sigma) field is gauged according to the prescription 
used in Refs. [g |. 

In the = gauge the static Hamiltonian is bounded from below by the sum 
of the two topological charge densities, the monopole charge and the degree of the 
map of the S 3 field on R3.TILUS the imposition of spherical symmetry reduces the 
system to an one dimensional subsystem, and the resulting differential equations are 
integrated analytically in the asymptotic regions and then numerically. This yielded 
the monopole-Skyrmion. 

In the Aq 7^ gauge, the Euler-Lagrange equations arising from the variation 
of the static Hamiltonian density do not yield a soliton with non-vanishing Aq and 
hence have Eq = 0. Instead, the variational equations arising from the (non positive- 
definite) Lagrangian density in the static limit support spherically symmetric solu- 
tions with Ef 7^ 0. This is also what happens with the JZ dyon. There 0], inspite 
of the non-positive-definiteness of the functional subjected to the variational princi- 
ple, it happens that after taking the static limit and imposing spherical symmetry, 
these equations reduce to a set of consistent, i. e. not overdetermined, set of cou- 
pled ordinary differential equations. Their solutions support a non-vanishing Aq 
field. These ordinary differential equations also result from the variation of a certain 
one-dimensional (radial) functional which, in contrast to the one dimensional energy 
functional, is not positiv definite. 

In the light of the surprisingly successful outcome for the JZ dyon, we were moti- 
vated to address the same question for the SO(3) gauged 0(4) model [§, [3[]. Subjecting 
the Lagrangian to the variational principle and then taking the static limit and im- 
posing spherical symmetry, we found that this also led to a consistent set of coupled 
ordinary differential equations. The same situation obtains with the composite model 
of this paper, and it is the dyon like soltions of these last equations which yielded 
the dyon-Skyrmion. Concerning the SO(3) gauged Skyrme model on its own, while 
its equations of motion reduce to a consistent set of coupled ordinary differential 
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equations, their solutions support only vanishing electric field. 

As a byproduct of our study of the dyon-Skyrmion, we made a detailed re-analysis 
of the JZ dyon refining our understanding of the latter, for example exploring the 
dependence of the energy of the dyon on its electric charge. 

An important result of the numerical analysis of the monopole-Skyrmion solu- 
tion is that, as the coupling strength of the Skyrme term is shrunk down to zero 
the monopole-Skyrmion reduces to the monopole, as depicted in Fig. 1. Thus the 
monopole does not stabilise the 5*0(3) gauged sigma model without a Skyrme term, 
something that is not prohibited by the Derrick scaling requirement. 

Perhaps the most interesting aspect of the dyon-Skyrmion occurs for the model in 
the PS limit (A = 0) in the special case where the boundary value q of the function 
g(r) parametrising Aq equals 1. In this case, the equations governing the functions 
a(r) (parametrising the gauge field) and the function f(r) (governing the Skyrme 
field) decouple from the fields h(r) (parametrising the Higgs field) and g(r). As a 
consequence the solutions for the functions a(r) and f(r) are just the (branch B 
of the) gauged-Skyrmion solutions and exist only for values of the Skyrme coupling 
constant larger than a critical value k%, as seen from Fig. 2. However, when the 
integrations of the Higgs field function h(r) and of the dyon function g(r) are taken 
into account, then finite energy solution only exist if the Skyrme coupling constant 
is larger than the critical value Ky 2 > k^, see Fig. 9. The energy of the solution 
at this critical value is found to become infinite and for lower values of the Skyrme 
coupling constant no finite energy solution exists. The time rate of change of the 
chiral charge eq. ( f44[) is equal to the integer 1 (up to normalisation) for all values 
of the Skyrme coupling constant k down to the critical value k^ 2 , below which no 
finite energy solutions exist. We hope that this result may prove relevant to the 
semiclassical description of monopole catalysis of Fermion number non- conservation. 
If for example it could be argued that the dyon-Skyrmion favoured by Nature is the 
solution to the system eq. (0) in the PS limit, with the asymptotic constant q — 1, 

i. e. for which the quantity is an integer (up to normalisation), then it would 

follow that below the critical value Ky 2 there will be no Q5 violating rate. We intend 
to return to this question in the near future. 
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A 


monopole 


monop ole- Skyrmion 
(£ = 1, « = 0.4) 


gauged Skyrmion 
(£ = !,« = 0.4) 


0.0 


1.000 


3.450 


2.98 


0.05 


1.106 


3.470 


2.98 


0.10 


1.138 


3.480 


2.98 


0.20 


1.180 


3.490 


2.98 


0.40 


1.220 


3.510 


2.98 


0.60 


1.250 


3.520 


2.98 


0.80 


1.270 


3.530 


2.98 


1.00 


1.290 


3.536 


2.98 



Table 1 

The energies of the monopole, the monopole-Skyrmion and the gauged Skyrmion for several 
values of the Higgs coupling constant A. 
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Figure Captions 



Figure 1 The energies eq. (25) of the monopole (line a), of the gauged Skyrmion (line b) and 



of the monopole-Skyrmion (line c) as functions of k (A = 1,£ = 1). 

Figure 2 The energy of the gauged Skyrmion as a function of k in the region of the phase 
transition. The branches A, A are represented by the solid line and branch B by the 
dashed line. 

Figure 3 The (logarithm of the) function a{x) on the two branches B and A on a logarithmic 
scale for several values of k approaching the critical value k c £. 

Figure 4 The quantity F\ defined in eq. (^) is plotted as a function of k for the branches A, A 
(solid line) and for the branch B (dashed line). 



Figure 5 The quantities ln(F — 1) and ln(a) (defined in eq. (p4|)) are plotted as functions of 
the parameter ln(«; — Kg). 

Figure 6 The values of the electric charge Q as a function of the parameter q. The solid lines 
represent the dyon for A = (line a) and A = 0.5 (line b). The dashed lines represent 
the dyon-Skyrmion (£ = 1) for A = 0, k = 0.4 (line c), A = 0.5, k = 0.4 (line d) and 
A = 0, k = 1 (line e). 

Figure 7 The values of the energy E as a function of the parameter Q. The solid lines represent 
the dyon for A = (line a) and A = 0.5 (line b). The dashed lines represent the dyon- 
Skyrmion (f = 1) for A = 0, k = 0.4 (line c), A = 0.5, k = 0.4 (line d) and A = 0, 
k = 1 (line e). The stars depict the points where the solution has maximal finite 
charge. The bullets correspond to the data of @. 

Figure 8 The value of the critical charge Q cr as a function of A. The solid line refers to the 
dyon solution. The dashed line refers to the dyon-Skyrmion solution for £ = 1 and 
K = 0.4. 

Figure 9 The energies of the dyon-Skyrmions with q = 1 (solid line) and q = 0.5 (dashed line) 
as functions of k (A =0, £=1). The vertical dotted line indicates the critical value of 

K. 
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Fig. 5 
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Fig. 6 
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